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1. Introduction 

Resonances in Quantum Mechanics describe states evolving non-unitarily [TJ, both 
decaying and building ones, and have found applications in the study of nuclear 
reactions. As such, they have been more often studied in the case of spherical nuclear 
potentials, mainly of the "shell-model" type and restricting the analysis to the radial 
s-wave equation. Albeit for obvious unessential differences in the Boundary Conditions, 
the physical picture is similar to the case of the 1-dimensional barriers we are considering, 
more akin to condensed matter systems or ion-trapping devices. The interest is alive 
and work is still in progress on the physical interpretation of resonances [2]. 

Resonant solutions to the Schrodinger equation occur in any simple potential 
barrier with a compact support, the plain square barrier being the most simple, 
tractable and fully representative example, and they are more often found in alternative 
specialized expansions of the Green Function and other physically relevant objects. They 
correspond to complex energies and momenta, which place them out of the familiar realm 
of Hermitian operators in Hilbert spaces, where a wealth of well known mathematical 
properties is available. Though having been first studied in early times [3]j4], the non- 
trivial mathematical properties of resonances have spurred a long lasting investigation 
effort. This paper focuses on some still debatable issues regarding the norm, inner 
products and completeness properties of resonant states. 

The issue of the completeness arises in its simplest form as soon as one attempts, for 
instance, to expand the identity in terms of projectors on bound (if any) \(j>i) , resonance 
\z n ) and background ("scattering" complex energy) \z) states, namely 

I = Y,\^^+H\ Z n)^n\ + hz\z)( Z \ . (1) 
i n J 

where the second sum may involve a variable number and type of resonances. Such an 
expansion is, in principle, attainable by deforming, in the complex (two-sheet) plane [5], 
the continuum real energy integration occurring in the traditional expansion in terms 
of bound and scattering states 

J = EI&>M + J dE \ E )( E \ (2) 

i 

following the lines of the proof of (2), as given for instance by [6]. In particular cases 
a complete expansion in terms of resonances (plus bound states) can be found, as for 
example for continuum wave functions (or the scattering solutions) within a finite region 
including the support of the barrier [7] [8] [9]. This has direct application for instance 
in time- dependent problems, as in the time evolution of quantum decay [10]. In the 
general case, covering the whole space, directly testing any resonance expansion like 
(1) against the idempotence requirement I 2 = I, needs the computation of the inner 
products of all these families of states with each other. We aim to perform some steps 
in this direction. 

The difficulty stems from the divergent asymptotic behaviour of the resonance 
solutions, which leads to infinite norms and seemingly divergent and hard to calculate 
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inner products, and has given rise to a variety of proposals to circumvent it. These 
proposals adopt different prescriptions to render finite the space integrals involved in 
these products and in general matrix elements. We may quote analytical continuations 
of the resonant solutions in the complex momentum plane [11], the "External Complex 
Scaling" [12] of the space-coordinate integration variable, and the introduction of 
convergence factors in the integrals. 

In this paper we shall adopt a Gaussian convergence factor, first introduced by 
Zel'dovich [13] and used by others [5], with which a limit can be worked out yielding 
novel non-trivial results: in fact we recover the result of the integrals that were already 
well defined, and extend the finite result to a newer region of the complex momentum 
plane. The procedure, which avoids relying on analytical continuation arguments, leads 
to a prescription for the definition of the inner products, yielding a specific set of 
orthogonality relations. Our results are derived for a general potential with finite support 
of which the bound and resonant poles of the S-matrix are known. 

The usual inner product in the Hilbert space, involving the complex conjugate 
of one of the wave functions, is defined so as to have real probability densities and 
norms for the general wave functions. However when dealing with resonant solutions, 
a symmetrical inner product is involved in the convolution of the Green Function with 
initial states (superpositions of scattering states), and is associated to a complex "norm" 
p3j, so we explore the properties of both alternative definitions. 

The plan of the paper is the following: In Section 2 we briefly review the relevant 
features of the resonances while fixing some notation and stressing the relationships 
between proper, anti, outgoing and incoming resonant solutions. In Section 3 we 
calculate the standard inner product of resonant states between themselves and with 
the bound and the scattering states. In Section 4 we introduce the symmetrical inner 
product, and give the result for the resonant-resonant and resonant-scattering product. 
Then in Section 5 we outline in some detail the regularization prescription adopted, 
based on a Gaussian factor, and show how the (either finite or divergent) results are 
recovered in the non- regularized limit. For comparison, we comment here on some other 
prescriptions found in the literature. Finally the Conclusions are drawn in Section 6. 
The crucial integral formulas are derived in an Appendix. 

2. The resonant solutions 

We consider the 1-D time-independent Schrodinger equation 

[— + p 2 - 2mV{x)\ i/t(x, p) = (3) 

where we use units such that h = 1. 

For a cut-off potential V(x) describing a general barrier with support in the compact 
interval [0,L], besides the usual scattering in and out solutions for continuous real energy 
E = p 2 /2m > and possible bound states (discrete real E{ < 0), one has resonant 
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solutions (Gamow states) satisfying the homogeneous Outgoing Resonant Boundary 
Conditions (ORBCs), 

d x ip\ x = = -ipV(O) » d x ijj\ x=L = vp'il)(L) . (4) 

Solutions w n (x), called proper (Outgoing) resonances, exist for a denumerable set of 
isolated values p n of p (with corresponding energies z n = p 2 l /2m) lying inside the 
octants close to the real axis (i.e. |Rep n | > |Imp n |) in the lower half complex plane 
(i.e. Imp n < 0), and occupy symmetrical positions with respect to the imaginary axis 
(see Figure 1). It is customary to label them as p n (n — 1, 2, ...) when Rep n > and 
as p_ n = — p* n their symmetric ones, some times called anti-resonances. In the case of 
a simple square barrier, the real parts Rep n tend to be spaced regularly for increasing 
\n\, while | Imp n | grows slowly [T]. 

RBCs with reversed sign of p in (4) correspond to Incoming solutions u n (x), and 
the momenta p n = — p n = p*_ n lie in the upper half complex plane. We denote with 
u n (x) = u n (p n ; x) these solutions and with \z n ) the corresponding states. Notice that the 
corresponding energies z n lie in the first Riemann sheet. For real potentials V(x), the 
complex conjugates u* n (x) of the outgoing resonances correspond instead to yet outgoing 
solutions with mirror momenta — p* , and energies z* n in the second Riemann sheet, so 
that u* n (x) = U- n (x) and z* n = Z- n . For the same reason, the complex conjugates of the 
incoming resonances are again incoming solutions. In the literature it has been easily 
mistaken u* n (x) for an incoming solution. 

The ORBCs are equivalent to imposing the asymptotic form 

. / R n e~^ x , x<0 ... 
Un{x) = \T n e^ , x>L (5) 

where the amplitudes R n and T n differ by a phase and are defined up to a global arbitrary 
factor. An immediate consequence is that the norm 

/+oo f+oo 
di M*(j;)« n (x) ~ / dxexp2|x||Imp n | = oo (6) 
-oo J —oo 

is even more divergent than for the Dirac states ||^|| = (E\E) = 6(0) ~ f_™dx = oo, 
causing both kinds of states not to belong to C 2 . 

This fact, as opposite to the finite norm of the bound states, has been addressed in 
the framework of Rigged Hilbert Spaces (RHS) pJ)][T6] along the same lines adopted for 
the Dirac kets. In RHS, these states are interpreted just as linear functionals on more 
restricted spaces (Hilbert or Schwartz) in which an inner product is properly defined, but 
nevertheless it is customary to consider, in a somewhat relaxed sense, inner products of 
the Dirac kets which have a meaning only as distributions, namely (E\E') = 5(E — E'). 
However, when one tries to give an answer of this type for the resonance states, the 
inner products 

f+OO 

(Zn\Zm) = dx U* n (x)u m (x) (7) 



oo 
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may be expected to be generally divergent as well, because of the exponential growth for 
large x . For this reason their actual calculation has spurred the adoption of a number 
of elaborated strategies. 

The same difficulties affect the inner products (z n \E) between resonant and 
scattering states and, more generally, between resonant states and wave packets. This 
is an important issue as long as the analytical calculation of, for instance, the time 
evolution of an initial state impinging on a barrier [9] , or of the shutter problem [17] , 
contains sums of terms with poles in the resonant momenta. For the (either Laplace or 
Fourier) transform of the time dependent Green Function one has [1] 



G(x, x'\ p) = T-r t; ^ n ^ — - + regular terms (8) 

V N n 2p n p-p n 

The interpretation of the pole terms as (transient) excitations of the resonant eigenmodes 
of the system requires the knowledge of the projection of the initial state on the resonant 
ones. The dimensionless normalization factor in (8) is |14]|17] 



where a = and b = L, which is a sort of (complex) "norm". Notice that, when applied 
to wave functions ?p(x) £ C 2 , it reminds the usual (real) norm in Hilbert space albeit 
for replacing the squares by the square modulus and letting a and b recede respectively 
to — oo and to +oo. For the resonant solutions however, N n diverges in this limit. 



3. Inner products 

In this section we shall calculate the usual inner products (z n \z m ) of resonant states 
with themselves (7) and their mixed products (z n \E) and (<f>i\z n ) with the scattering 
and bound states respectively. The method leads also to the usual inner products 
= , (<t>i\E) = , and (E\E') = 5(E - E'). 



3.1. Resonances 

The space integral giving the inner product of two resonant solutions can be split in 
three sectors 

(z n \z m ) = dx U* n (x)u rn (x) 

= R* n R m dx e ^-Pm> + jL dx u * nUm + T * Tm dx e -KPn-Pm)x ( 10 ) 
= {R* n R m + T*T m ) / °° dx e-Wn-Pm)* + jL dx u>m _ T*T m J L dx e^-p™)* 

The ORBCs in and L, together with their complex conjugates, allow to express 
the finite integral J L dx u* n u m in terms of the amplitudes R and T outside the barrier, 
regardless of the explicit form of u n (x) inside the barrier and hence of the particular form 
of the potential. The procedure uses the Schrodinger operator O = —2mH = d 2 — 2mV 
and integration by parts to obtain 

(P*n-Pm)f dxu* n u m = 2m{z* n - z m ) [ dxu* n u m (11) 
Jo Jo 
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f L dx u* n (x)0 - t)u m (x) = W[u* n , u m ]% 
Jo 



where W[(p = <p d x tp — ipd x </> is the Wronskian of the functions <fi and ip , so that 



dx u* n u m 



jrp*rp -i(p*-p m )L , D* D 



(12) 



'0 K - Pm 

Then cancelations occur such that 

POO i 

(z n \z m ) = (R* n R m + T*T m )( / dx e-^-nJ* + ) (13) 

Now the point is the calculation of the integral in (13). We adopt the limit A — > 
in (29) as a prescription for the result of (13), and defer the discussion of our scheme to 
Section 5 and to the Appendix. 

With this prescription we finally obtain 

'0 , -f < arg(p m -K) < 5f 
oo , otherwise 



(14) 



In particular (z n \z n ) = oo, as expected. 

For each resonant state \z n ), the result (14) defines a "neighborhood of divergence" 
so that \z n ) is orthogonal to any other \z m ) the momentum p m of which lies outside a 
"cone of divergence" with apex in p* , and gives a divergent inner product if p m lies 
inside this cone (Figure 1) 
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Figure 1. An example of orthogonality between resonant states: 
(z2\z n ) = for n < and n > 4 , while (z2\z m ) = oo for < m < 4 . 
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3.2. Resonant and Scattering states 

The in and out scattering solutions obey a single differential BC at one of the points 
x = or x = L . For instance, a right-moving in state obeys only the second BC in (4) 
with p > 0. This is equivalent to imposing the asymptotic form 

^i^e- , x>L (15) 

For instance, the combined BCs of the resonant and of the (p > , m) scattering 
solutions similarly lead to 

POO 

(z n \E)= / dx{R* n e- iip « +p)x + (R* n R + T:T)e- iip "- p)x ) (16) 
Jo 

-(i^ + T n *T) + — —R* n 



Pn - P Pn+P 



With p > , for n > we always have — | < arg(p* + p) < 5| , so that 

(^|£> = (i^ + ^TX rdxe-**-& + V— ) (17) 



(*n|£> 



and therefore, with the prescription adopted, 

"0 , -f < arg(p-p;) < 5f 
oo , otherwise 

This means that a given scattering in state (with momentum p > on the real 
axis) is orthogonal to any \z n ) (with n > 0) if the momentum p n lies outside the cone 
with apex in p , the inner product being divergent otherwise. Viceversa, given p n , 
the momenta p of the orthogonal scattering states lie outside the cone with apex in p n 
(Figure 2). 

The scattering states always have a reflected wave with momentum of opposite sign 
to the incident one, so the situation is trickier for n < . Given p-\ n \ , the cone to be 
considered is again the one with apex in the mirror momentum p\ n \ . 

3.3. Resonant and bound states 

Let us suppose that, besides the barrier, there is some potential well within the region 
[0, L] sustaining bound states \4>i) with purely imaginary momenta pi = iqi (g, > 0). The 
exponential decrease of the amplitude outside the well, corresponding to resonant-like 
BCs 

d x (p\ x=0 = q 0(0) , d x (t)\ x=L = -q (f>(L) (q > 0) , (19) 

which imply the more popular (and weaker) 0j(±oo) = 0, may be translated into the 
assumption of an asymptotic form similar to (5), namely 

Hx) = \ 1 ' (20) 

r v ; \ Tie-** , x > L K J 

which is manifestly square integrable. 
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Figure 2. Two examples of orthogonality between resonant and scattering states: 
(z2\E) = (oo) when p = \/ 2mE C R+ lies out of (inside) the interval [A2,B2\ , 
and {E c \z n ) = (oo) when n ^ 3,4 {n = 3,4) respectively. 



Following the same steps leading to (13) we now have 

r°° i 

(<j> i \z n ) = (R*R n + T*T n )( dxe^^-— ) (21) 

Jo iqi + p n 

The integral is convergent, yielding i(igj +p n ) _1 and hence (<fii\z n ) = , provided that 
qi > | Imp n | . Thus we have a situation similar to that occurring between the resonances, 
namely that the states are orthogonal if the bound state momentum, lying in the positive 
imaginary axis, and the resonant one, lie outside the respective divergence cones, the 
inner product being infinite otherwise. It might happen that, for general analyticity 
reasons for any general potential V(x) with barriers and wells within [0, L], the cones 
of the allowed bound state momenta do never include the resonant momenta, but this 
point would require separate investigation. 

4. Symmetrical inner products 

Together with the standard definition (7), which yields (finite or infinite) real norms, 
we shall consider also the alternative symmetrical definition 

/+oo 
dx <p{x)ip{x) , (22) 
-oo 

which arises in the convolution of the kernel (8) with the initial state, and yields complex 
"norms" . 
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{z n \z m } = (R n R m + T n T m )( / dxe^ +p ™> ■ ) (23) 

JO V n + Pm 



The result is 



w«.} = {° • :t<^(p m +p.)<5f 

oo , otherwise 



For a given \z n ), the apex of the divergence cone is at — p n so that, in particular, 
{z n \z m } = if both m, n > 0, also if both m, n < 0, and, more noteworthy, {^nl^} = 
for any n . Only some of the products are divergent when n and m have opposite sign. 
For the crossed products of \z n ) (n > 0) with the scattering p > in state one has 

{z n \E}= / dx^Rne^-^ + {R n R + T n T)e^ + ^ x ) (25) 
Jo 

-(R n R + T n T) _R n 



Pn+P Pn-P 



One always has — f < arg(p n + p) < 5^ , so that 

K|£} = { ° ' ~ ^ < ar s(P« -P) < 5 I (26) 
n [ oo , otherwise 

Therefore the divergence cone has the apex in p n , as in the case of the standard inner 
product (z n \E). For n < , only p + p_|„| may lie in the sector 5^ < arg(p n +p) < 7j 
leading to a divergent product, so that, for a given p~\ n \ the p -states are orthogonal if 
p lies inside the cone with apex in the mirror momentum p\ n \, reproducing the same 
situation encountered for the standard product. 

Because of the relationship U- n (x) = u* n (x) between the outgoing solutions, we 
have (z n \z m ) = {z_ n \z m } and {z n \z m } = (z_ n \z m ). Then the result {z n \z n } = is 
less surprising if rewrite it ) = , just a particular case of the conventional 

orthogonality between states with momenta in opposite quadrants of the lower half 
complex plane. 



5. Regularization of the divergent products 

Several proposals have been worked out in the literature to deal with the divergent 
inner products. We quote here the one by W. Romo [11] of analytically continuing the 
momentum dependence of the solutions u(x) inside the products from the upper-half 
complex plane, where the (outgoing) function would be square-integrable, to the lower 
half-plane where the resonant momenta lie. This is equivalent to prescribe the finite 
value ik^ 1 for the result of the integral Jq 00 dx e lkx also in the whole lower half-plane k, 
where it would actually be divergent (see the discussion in Appendix). This approach, 
combined with advantages of working in momentum representation, has been followed 
in later works |18j . 
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Other proposals rely on the introduction of convergence factors in the integral which 
are able to cope with the asymptotic exponential growth of the resonant solutions. The 
simple factor e~ Xx , (A real > 0) works only if A is greater than the absolute values of the 
imaginary parts of given p n and p m (which grow with increasing n , m), so that A cannot 
be brought to zero without rendering divergent the integral. However it lets performing 
an analytic continuation in A to its "forbidden" values which may be adopted as a 
prescription. 

Following Zel'dovich [13] , we adopt the more powerful Gaussian factor e~ Xx2 , which 
is able to overcome the growth of all the resonant solutions even for A — > + , and is 
computationally tractable |13j[5]: 

r+°° \ 2 

({z n \z m )) x = / dx e~ Xx u* n {x)u m (x) (A real > 0) (27) 

J — oo 

The functions e~^ Xx2 u n {x) are not eigenfunctions of the Hamiltonian so the exact 
calculation of (27) requires the knowledge of the explicit form of the solutions inside the 
barrier. 

Equation (10) now reads 

POO 

((z n \z m ))x = (KRm + T* n T m ) / dx e^V^"^ (28) 

Jo 

+ [ dxe- Xx \* n u m -T*T m f dx e~ Xx2 e- i{p -- p ™ )x 
Jo Jo 

According to (A.l), this regularized expression is finite for A > : 

((z n \z m ))x = (KR m + KT m ) ^ V^ze z2 erfc(z) (29) 

Pn Pm 

+ ( L &xu* n u m -T* n T m [ L dxe- iip «- p ™ )x +0(A) 
Jo Jo 

where z = \(p* n - p m )/ (2y/X) . 

In the limit A — > + the terms 0(A) vanish, and for y/ixze z erfc(z) we have the 
result (A. 5). Thus we obtain (14) in this limit, which must be intended as a prescription 
for the calculation of (13). We discuss this issue in more detail in the Appendix. 

Notice that for finite A > 0, the exact integrals over the finite interval [0,L] get 
rather involved and the one over the resonant functions inside the barrier depends on 
the explicit form of the functions other than on the position of the resonant momenta. 

The regularization works in a similar fashion for the crossed inner products (z n \E) 
and (z n \(f)i), as well as for the symmetrical definition. 

6. Conclusions 

We have calculated some relevant inner products involving the resonant eigenstates 
in the example of a 1-dimensional quantum potential barrier with compact support. 
This old problem is non-trivial since the modulus of the resonant solutions grows 
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exponentially at the spatial infinity, giving rise to infinite norms and seemingly infinite 
inner products. 

Among the variety of historical proposals to circumvent these difficulties, we have 
adopted a Gaussian convergence factor, first introduced by Zel'dovich, and carried out 
the limit of the integrals where the factor fades off to unity This prescription yields 
inner products such that most of these states are orthogonal to each other, except when 
they lie in a neighborhood defined by a "divergence cone", in which case the product 
is infinite. Similarly, scattering states (with real momentum p > 0) and bound states 
(pure imaginary momenta iq^ , with ^ > 0) are orthogonal to a given resonant state 
\z n ) with momentum p n , except when p or iqi lie in a neighborhood of p n . Thus, 
the resonant states share properties of the continuum and discrete real spectra, namely 
(partial) orthogonality and the infinite norm characteristic to the Dirac states. 

This result is different from the full bi-orthogonality obtained by the prescription 
of analytically continuing the finite integrals from the upper complex momentum plane 
to the whole lower half plane, where the resonant momenta lie but where the integrals 
are formally divergent. Our limiting procedure instead extends the finite result to the 
more modest 7r/4 angular sectors of the lower half plane close to the real axis. We differ 
also from earlier attempts [5] using the Gaussian convergence factor, actually limited to 
the products (z_ n \z m ) (where \z- n ) was mistaken for the incoming state \z n )), in that 
our exact calculation yields also (z_ n |z n ) = instead of a finite quantity normalizable 
to 1. It is interesting also the comparison with the results obtained in the Friedrichs 
model PS] [2D]. 

Therefore the orthogonality and normalization properties obtained depend on the 
prescription adopted, although one should expect, as a signature of their consistency, 
that all of them lead to the same unique result. For instance, the inner self-product 
(z n \z n ) yielding the square of the conventional norm \\z n \\, which is manifestly divergent 
on the same footing of the scattering states, is finite for some prescriptions, whereas 
ours recovers the infinite result. 

We do not investigate here the question if a regularization by different convergence 
factors (using for instance e~ x ^ for some real value 1 < v < 2) would yield narrower 
divergence cones, thus approaching the result (z n \z m ) = (n ^ m) which would seem 
more natural and closer to the result of the analytical continuation prescription. On 
the other hand, such a continuation in the momentum of the solutions is a non-trivial 
matter [21] related to the two- sheet structure of the complex energy plane. 

The completeness of the scattering (plus bound, if any) solutions in Hilbert space 
is not inherited by the sole resonant (plus bound) states in resonance expansions of 
the unity, and a continuum of complex energy "background" states must be included 
[5] [22] [23]. Any consistent unambiguous prescription for the calculation of the inner 
products involving resonances should be tested against the requirement of idempotence 
of the unity, as discussed in the Introduction. It is not clear how could it happen 
in the analytical continuation prescriptions, where the claimed orthogonality between 
resonances (or between resonant and anti-resonant states) is of the Kronecker-5 type, 
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and the cross inner products between bound, resonant and background states have not 
been worked out. 

We instead have orthogonality between most of the states and divergent inner 
products between "neighbor" resonant (and continuum and bound) states, which then 
have also infinite self-products consistently with the well known infinite norm of the 
resonances. In a general expansion involving bound, resonant (in variable number, 
according to the poles encircled by the deformed integration path), scattering and 
continuum background states, showing directly the idempotence of the unity requires 
tackling the products involving all of these families of states. In this paper we have done 
part of this program and developed some regularization and calculational techniques. 
Work is in progress on the inner products involving the resonant and the background 
states, aiming to the eventual cancelation of the divergences encountered. 
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Appendix A. Limits of erfc( z) -related integrals 

We rely on the basic integral 

J(fc, A) = / dx e~ Xx e lkx = - ^ze z eric(z) (A real > 0) (A.l) 
Jo k 

which is directly related to (7.1.2) in [21], where z = —\k/(2y/X) , and hence k , can take 
any complex value. Notice that on completing a square in the exponent and changing 
the integration variable to t — \/\ x + z one obtains the intermediate expression 

1 „ fOO+Z „ 

J(k,\) = - r e z die-* (A.2) 

V A Jz 

The integral representation (7.1.2) is convergent when the path of the complex 
integration variable t approaches oo along a direction — f < arg(t) < | . This condition 
is fulfilled by (A.2) as long as z remains finite. However, in (A.l) not always the limit 
A — > (and hence z — > oo) can be taken before the integration is carried out, because 
the function f\(x) = e~ Xx2 does not converge uniformly to the function fo{x) = 1 when 
A . 

For Im k > the limit can be taken in the integrand in (A.l) because the finiteness 
of the integral is always assured by the convergence factor e ~( lrnk ) x ^ trivially yielding 
i k~ x . In this case J(k, A) converges to J(k, 0) = ik~ l when A — > . 

For real k the result 

i 

dxe ikx = iPV- + 7r5(k) (A.3) 
k 
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is also well known. Notice that it relies on adding to k a small imaginary part ie, which 
still guarantees the convergence when A — > 0, but later in the limit e — > + the result 
must be interpreted distribution. 

For Im k < the integration and the limit A — > do not commute. In that case 
we adopt the limit of the integral as a prescription. In the right-hand side of (A.l), the 
limit A —7- can be directly inferred from the asymptotic expansion 

I- 7 2 r , \ / \m 1 ' 3 ••• (2m - 1) 37T . . . 

v^e 2 erfc( 2 )~l+^(-l) m ^ (z^oo, |arg 2 | < — ) (A.4) 

m=l v Z2; / 4 

(see (7.1.23) in [21]), namely 

r- z 2 r r s f 1 , ~3f < arg(z) < 3f , . r , 

hm^ooA/TT^e 2 erfc(^) = { , 4 w 4 (A.5) 

[ oo , otherwise 

which can also be numerically checked. In this way (A.5) yields an extension of J(k, 0) 
to a new region of the lower half complex plane: 

J(*,0) = f* ' i <arg « <5 5 ■ k *° (A.6) 
oo , otherwise 
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